Rhythmic cluster generation in strongly driven colloidal dispersions 



H. H. Wensink and H. Lowen 
Institut fur Theoretische Physik, Heinrich-Heine-Universitat Diisseldorf, 
Universitatsstrafte 1, D-40225, Diisseldorf, Germany 
(Dated: February 6, 2008) 

We study the response of a nematic colloidal dispersion of rods to a driven probe particle which 
is dragged with high speed through the dispersion perpendicular to the nematic director. In front 
of the dragged particle, clusters of rods are generated which rhythmically grow and dissolve by 
rotational motion. We find evidence for a mesoscopic cluster-cluster correlation length, independent 
of the imposed drag speed. Our results are based on non-equilibrium Brownian dynamics computer 
simulations and in line with a dynamical scaling theory. 
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In recent years particle tracking techniques have been 
developed to determine the viscoelasticity of soft mate- 
rials on a microscale Q. Optical or magnetic tweezers 
can be used to drag a bead particle through the mate- 
rial in a controlled way. The system's response to the 
tracer particle reveals its microrheology and provides in- 
sight into the viscoelastic response of biological matter 
U , colloidal dispersions 0, 0] , polymer solutions and 
wormlike micelles Q. 

Microrheological measurements are usually based on 
linear response, employing the intrinsic thermal motion 
of the probe particle. Much less is known about non- 
linear effects that occur when soft materials are subject 
to strongly driven probe colloids as expressed by a high 
Peclet number. The best studied examples are colloidal 
suspensions thanks to their model character and direct 
accessibility by real-space techniques 7]. Striking novel 
effects in strongly driven colloidal dispersions include 
force-thinning of the microviscosity ||, local jamming of 
particles in front of the dragged particle [jj] , and rhythmic 
bursting of colloids driven through channels |ic| . 

In this Letter we propose a new effect of colloid cluster 
generation by a spherical probe particle which is dragged 
through the suspension with a high velocity. In contrast 
to static colloid clustering which was recently observed 
under equilibrium conditions 0, i the cluster genera- 
tion is rhythmic here, i.e. particle clusters grow and dis- 
solve with a characteristic frequency. Correspondingly, 
along the trace of the moving probe particle, a character- 
istic cluster-cluster correlation length Ao occurs. While 
the average cluster size increases monotonically as a func- 
tion of the drag velocity, the correlation length is found 
to be independent of the velocity for strong drags. 

Our results are obtained by Brownian dynamics com- 
puter simulations of a two-dimensional (2D) suspension 
of rod-like colloidal solute particles. The rods are in the 
nematic phase and the probe particle is dragged perpen- 
dicular to the nematic director and cannot penetrate the 
rods. We observe the following cluster generation pro- 
cess at high drag speeds: the moving probe particle will 
touch a rod which is then taken up to its speed forming a 



particle-rod complex. This complex successively sweeps 
together further rods such that a cluster grows in front 
of the dragged particle. Forcing a collective rotation of 
the whole cluster, the tracer particle is rushing over the 
cluster, the cluster slowly dissolves and the process starts 
again. In order to quantify and predict the statistics of 
this process of dynamical cluster generation, we propose 
a simple dynamical scaling theory. The simulation results 
for the cluster size distribution, the time-averaged cluster 
size and the characteristic cluster correlation length Ao 
are in agreement with the theory. 

The rhythmic cluster generation effect is very general. 
It is also expected in three dimensions and for tracer and 
colloidal particles of different shapes. The only condition 
is that the interaction between the particles possesses an 
excluded volume region, i.e. a region of impenetrability. 
Therefore, rhythmic clustering should be observable in 
many different experimental set-ups. However, cluster 
generation is absent for a soft spherical tracer particle in 
a suspension of soft colloids as simulated in Ref. 

An important implication of the cluster generation oc- 
curs for suspensions which coagulate once their distance 
falls below a critical value. Since the average size of the 
clusters can be tuned by the drag speed of the sphere, 
it is possible to generate aggregates with tunable size. 
These may serve as building blocks for nano-composite 
materials with novel rheological or optical properties . 

In our Brownian dynamics simulations, we consider a 
2D system of N charged rodlike particles kept at constant 
temperature T via the solvent. The simulation involves 
a numerical solution of the overdamped Langevin equa- 
tions using the (short-time) rotational diffusion constant 
Dq and translational ones Dq and Dq , parallel and per- 
pendicular to the rod axis, for a single non-interacting 
rod with length Lq and thickness do [lfij . The mu- 
tual rod interaction is represented by a Yukawa-segment 
model in which the total rod charge is equally partitioned 
over n — 13 equidistant segments located on the rod 
axis. The interaction potential between two segments 
a and b on rods i and j (J ^ i) is of the Yukawa type, 
U = Uq cxp(— Kr^ b )/nr^ b with k the Debye screening con- 
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stant and r^ b the segment- segment distance. The am- 
plitude is U = 6.5 x 10~ 2 k B T (with k B T the thermal 
energy) comparable to that of Tobacco Mosaic Virus par- 
ticles in a salt-free solution under standard conditions for 
which Lq = 300 nm, cIq = 18 nm and kT x w 50 nm |l6|. 

A starting configuration was generated by randomly 
putting N = 900 rods in a parallel, non-overlapping con- 
figuration such that the nematic director is parallel to 
the y-axis of a rectangular simulation box of lengths L x 
and L y . In all simulations the reduced number density of 
the nematic system was fixed at p* = NLq/ L x L y = 3.0. 
Equilibration typically required 10 5 Brownian timesteps 
of magnitude At = 0.002rJ, measured in units of the 
typical translational Brownian relaxation time of a free 
rod r T = L§/r>g\ with D% = (Z)jJ + 2L» ( f)/3. We con- 
firmed that this system is in a two-dimensional nematic 
phase [l7l | with a director along the y-axis of the simula- 
tion box. 

The next step in our simulation is to drag a small 
sphere through the nematic host structure along the x- 
axis of the frame. The probe particle is a small repul- 
sive sphere interacting with the rod segments via a trun- 
cated and shifted Lennard- Jones potential (setting the 
Lennard-Jones parameters to <tl,t = « and £lj = k B T) 
and a cut-off distance 2 1 / 6 ctlj- The particle is trapped 
in a parabolic potential well which is displaced through 
the system along x at a constant speed v. The external 
potential for the sphere located at position r can be writ- 
ten as Z7 ex t(r,i) = Aq\t — uix| 2 with a large amplitude 
A = IO^bT/Lq for the trapping potential. The par- 
ticle is dragged until it has traversed the horizontal box 
size L x . Henceforth, we will use a reduced drag speed 
v* = vLo/Dq which expresses v in units of the Brown- 
ian speed of a freely diffusing rod. Quantitative results 
were collected by averaging over consecutive independent 
initial configurations for the nematic state. 

In Figs, la-c three simulation snapshots are depicted 
corresponding to a dragging event at increasing drag 
speeds. At high drag speeds a rhythmic clustering effect 
is observed, as reflected in Figs, lb and c, in which rods 
are pushed together into an anisometric cluster which 
subsequently dissolves by collective rotation under the 
influence of the force exerted by the dragged sphere. This 
process is depicted explicitly in a time series of snapshots 
shown in Figs. ld-f. In all snapshots, voids can be ob- 
served in the wake of the dragged particle which slowly 
dissolve over time. The length of the void trail grows con- 
tinuously with drag speed. The same holds for the void 
size indicating that the time needed for the system to 
relax back to equilibrium also increases with drag speed. 

The instantaneous cluster size can be monitored us- 
ing a simple cluster criterion based on the distance of 
closest approach between rod pairs located in front of 
the drag sphere. A rod pair is assigned to a cluster 
if this distance is smaller than 0.5 times the average 
distance of closest approach in the equilibrium nematic 




FIG. 1: Snapshots of the nematic structure at (a) v* = 
100, (b) v* = 1000 and (c) v* = 4000. The driven sphere is 
indicated by the encircled dot. The sequence (d)-(f) depicts 
the growth and rotation of a single cluster at v* = 4000. 
The length of the arrow corresponds to the cluster-cluster 
correlation length of Ao = 5Lq. 
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FIG. 2: Time series for the cluster size N c for a single drag 
event. Inset: time averaged cluster size (N c (vt)) versus drag 
speed. The dashed line {N c ) — 17.2 follows from theory. 



phase. The rhythmic nature of the clustering is clearly 
reflected in the cluster time series depicted in Fig. 2. 
Intermittent behavior is evident from the sudden out- 
burst of a big clustering event after a sequence of smaller 
ones. We verified that N c (vt) does not significantly 
change upon small variations of the cluster selection cri- 
terion. Defining (..) as an average over time t we ob- 
serve that (N c (vt)) increases smoothly as a function of 
the drag speed implying that the transition towards a 
state of rhythmic clustering occurs continuously. To in- 
vestigate correlations between the cluster sizes we have 
carried out a spectral analysis of the cluster time series. 
The power spectrum S(q) represents the Fourier trans- 
form of the cluster-cluster autocorrelation function via 
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FIG. 3: (a, left) Power spectrum for v* = 1000. The cor- 
relation peak is indicated by the arrow, (b, right) Dom- 
inant cluster-cluster correlation length Ao corresponding to 
the maximum amplitude of the power spectrum versus speed. 
The horizontal line is the theoretical result Ao ~ H.6Lq. 



S(q) = JZo dt ' exp[-iqvt'](N c (vt)N c (v(t + <'))) with q 
the wavcnumber. In the relevant spectral range, all spec- 
tra show a dominant peak (as indicated in Fig. 3a) which 
we can identify with a typical cluster-cluster correlation 
length X . Whereas the dominant peak is blurred some- 
what by the high noise level in the spectra at low speeds, 
it becomes much more pronounced at v* > 1000, see Fig. 
3a. The location of the dominant peak for various v is 
presented in Fig. 3b. A striking effect is that Ao does 
not increase with speed but remains at constant value of 
roughly 5 — &Lq. This result suggests the existence of 
a typical cluster-cluster correlation length scale for the 
nonlinear response of a driven system of colloidal rods. 
The length scale is also indicated in Fig. lc. 

Given the generic character of the clustering mecha- 
nism, we expect this length scale to be insensitive to the 
details of rod pair interactions and to depend only on the 
concentration p*oi the nematic phase |l8|. Moreover, we 
anticipate dynamic clustering to be unique for rodlike 
systems and do not expect it to occur in dense systems 
of soft spheres. The reason for this is that two charged 
rods in 2D have a finite core excluded volume in the limit 
\ 0, namely V eK = 2iQsin7 (with 7 the angle be- 
tween the rods), unlike soft Yukawa spheres for which 
V cx = 0. A finite exclusion zone around each particle 
is a necessary condition for the formation of large force 
chains that are associated with the clustered particles. 
To verify this statement we have performed similar drag 
experiments for dense fluids of Yukawa spheres close to 
the glass transition for various Lennard-Jones parameters 
for the sphere-rod interactions. In all cases, no evidence 
for rhythmic clustering was found up to drag speeds ex- 
ceeding the present range by several orders of magnitude. 

We now propose a simple dynamical scaling theory for 
the observed behavior. At high drag speeds the external 
forces due to the dragged spheres are much larger than 
the Brownian forces and the system is governed by over- 
damped, deterministic motion. The rotational motion of 
a single rod-like cluster with length £ under an external 
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FIG. 4: Sketch of the cage model used in the scaling theory. 



torque T is then given by £ r (£)ip(t) = T with ip the time 
derivative of the angle between the rod axis and y and 
£ r = £ J-(p) the rotational friction coefficient of the clus- 
ter, with T a function of the cluster aspect-ratio p [l9j . 
The force F directed along the x-axis is proportional to 
Icostp and the velocity, so we write F = v£j_£cos<p, with 
£j_ the translational friction coefficient of an individual 
rod. The torque is then given by the force times the ver- 
tical component of the arm length d. We consider d to 
be a dynamical angle-dependent variable proportional to 
speed, i.e. d = v sin ip. Finally, we assume the growth of 
the cluster length to be proportional to the speed while 
an increased inhibition is expected as the cluster rotates 
towards the x-direction. This gives i — vL$a~ l cos (p, 
where the additional length scale a pertains to the nearest 
neighbor distance in the nematic phase. It can be quan- 
tified from a simple cell description as sketched in Fig. 
4a. If we assume each rod to be confined to a rectangular 
cage, we can directly obtain a* = a/ Lq from the concen- 
tration of the nematic phase, via 2a* = (1 + 4//3*) 1 / 2 — 1, 
giving a* — 0.264 for p* — 3. We also assume the cluster 
aspect-ratio p to be a constant inversely proportional to 
a* , i.e. p = 1/a*. Gathering all contributions gives a set 
of three coupled dynamical nonlinear equations: 



v^ip =Q{x)^ L ^ 1 r 2 dcostp 
V d = simp 



* — 1 
a cos <p 



(1) 



with initial values {(^oj ^0} a t t — 0. It is important to 
note that the velocity can be scaled out of Eq. Q using 
the total drag distance vt as the independent variable. 
A stepfunction 0(x) with x = (t/2) costp — d has been 
included for consistency reasons to enforce the rotation 
and growth process to stop once the arm length exceeds 
1/2. 

Due to the external force the rod is pushed toward the 
edges of the cell, indicated in Fig. 4b, and we may define 
a typical cage deflection angle given by ip m = sin -1 a* and 
associated arm component d*j = (a* + l— cosy? m )/2. If we 
use these as the initial values for Eq. (JJJ we numerically 
obtain a dominant cluster length scale of Ao = vt(ip = 
0) = 14.6-Lo in qualitative agreement with the simulation 
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FIG. 5: Distribution of the cluster size N c for large speeds. 

data, see Fig. 3b. The overestimation is probably due 
to the somewhat crude nature of the cage model and the 
neglect of Brownian fluctuations which are expected to 
play a role in the initial stages of cluster formation. 

Fluctuations around the dominant length scale can be 
introduced by noting that the rods are not all parallel due 
to local orientational fluctuations around the nematic di- 
rector y. Therefore, we can introduce an angle 9 C measur- 
ing the deflection angle of the cage from its parallel ori- 
entation 9 C — 0. The fluctuation strength can be roughly 
quantified by a nematic order parameter < S z < 1, de- 
fined as S z = cos 20i/N) with < 0j < tt/2 the angle 
between rod i and y. We measured S z rs 0.8. Assuming 
the cell orientation angle 9 C to obey a Gaussian distri- 
bution P(9 C ) oc exp[— #e/2(l — S z ) 2 ], we can construct a 
series of cluster events by solving Eq. 0) for TV = 1000 
consecutive sets of initial values {</?,„ + 9 c ,d' m (9 c ), L }, 
with < 9 C < tt/2 sampled randomly from the Gaus- 
sian distribution (see also Fig. 4c). The cluster size N c 
can be obtained by multiplying the area fraction of the 
cluster rectangle located ahead of the drag particle with 
the typical rod concentration p* inside the cluster, which 
we estimate at p* « 6.0 from our simulations. From the 
time series (included in Fig. 2) we can extract the average 
cluster size (N c ) (Fig. 2) and the distribution function of 
the cluster size shown in Fig. 5. It can be deduced from 
Fig. 5 that the strong tail of the distributions at large 
speeds is correctly reproduced by the theory. 

In summary, we have studied the effect of a small col- 
loidal sphere driven through a nematic background at 
various drag speeds. We observe the formation of rod 
clusters whose growth is inhibited by collective rotation 
of the cluster under the influence of the force exerted by 
the dragged sphere. The growth-inhibition mechanism 
gives rise to rhythmic clustering behavior, quantified by 
a typical cluster-cluster correlation length which is inde- 
pendent of the imposed drag speed. A simple dynamical 
model without adjustable parameters is proposed based 
on overdamped rotational motion of a single rod clus- 
ter. The theory predicts a velocity- independent scaling 



value for the cluster-cluster correlation length in qualita- 
tive agreement with the simulation results. 

The observations described in this Letter can be ver- 
ified in experiments. Various experimental realizations 
are conceivable: one can confine suspensions of rod-like 
particles between two parallel glass plates and use an 
optical tweezer to drag a probe particle through the sus- 
pensions. Confocal microscopy and video-microscopy can 
be used to follow cluster generation dynamically. Fur- 
ther set-ups to test our predictions are colloidal nematic 
monolayers on substrates |20|. nematics anchored at in- 
terfaces |ilj, quasi-2D granular rods lip or elongated 
dust rods levitated in a plasma sheath 23] . Qualitatively 
similar rhythmic clustering is expected in 3D colloidal ne- 
matic systems if a thin transversely oriented rod is driven 
perpendicular to the nematic director. 
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